CO 



Dirac factorization and fractional calculus 

D. BabuscJ] 

INFN - Laboratori Nazionali di Frascati, 
Via E. Fermi, 40, 00044 Frascati (Roma), Italy 

G. Dattolfl and M. Quattrominfl 

ENEA - Centro Ricerche Frascati, Via E. Fermi, 45, 00044 Frascati (Roma), Italy 

P. E. Riccfl 

International Telematic University UniNettuno, 
Corso V. Emanuele II, 39 00186 Roma, Italy 



^ Abstract 

We show that the Dirac factorization method can be successfully employed to treat problems 
involving operators raised to a fractional power. The technique we adopt is based on an extension of 
the Pauli matrices and the properties of the roots of unity. We also comment about the possibility of 
using the method to linearize evolution equations containing the n-th root of differential operators 
and make a comparison with other techniques involving suitable transforms. 
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The formulation of the relativistic Dirac equation 111 is based on a factorization of the 

n 

square root that, as a consequence of the anti-commuting nature of Pauli matrices [2J 

{a j ,a k } = 2 6 jk (j,k = 1,2,3), (1) 

allows to write the Pythagorean theorem in the form 

VA* + B* = A&j + B& k (j^k), (2) 

where A, B can be either numbers or commuting operators. Eq. ([2]) can be extended to the 
case of three terms, and one can write 



VA 2 + B 2 + C 2 = Aa j + Ba k + Cai (j ^ k ^ I), (3) 



The realization of such a factorization is not unique. The Hamilton quaternions |3j can be 
efficiently used for this purpose, and this is by no means surprising if we note that a purely 
imaginary quaternion 

Q = ai + bj + ck, i 2 = j 2 = k 2 = ij k = -1, (4) 

can be realized in terms of Pauli matrices since one can write 

i = i&i, j = ia 2 , k = ia 3 . (5) 

It is therefore evident that the key tool of the game is the use of the Clifford algebras, at 
least if we limit ourselves to the case of square roots. 

We believe that a further element of interest in the use of this factorization method, and 
in its extension to higher-order roots, is that it can open new perspectives in the theory of 
the fractional calculus i. e., in the body of computations aimed at allowing the handling 
of operators (differential or not) raised to a non-integer number. Just to give an example, 
we note that, by choosing (j,k,l) = (1, 2, 3), we can write 



d z d x - % d y 



(6) 



d x + id y -d z 
while the operator Oi/ 2 = V d x + a can be cast in the form 

0i/2 = \fd~xbj + Vaa k (7) 



i.e., using the well-known identity |4] 
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For example, in the case j = 3, k = 2 one has 

1 



Ol/2 

or, by choosing j — 1, k — 2, 

1/2 = - 



—i a/7t ax 
71 x \ i y/irax — 1 



(8) 
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In a previous note {(3| we have shown an extension of the Dirac factorization method 
allows the extension of the identity ([2]) to cubic roots, namely 



vM 3 + B 3 = A Ai + B A 2 



where 
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(12) 



(13) 



with efc = e t27Tk / n (k = 0, 1, 2) cubic roots of unity satisfying the conditions 



Je fc = eo = 1, 



(14) 



fc=0 



fc=0 



from which we can show the following identities, necessary for the fulfillment of Eq. (TT21 
hold 

A| = 1, X 2 j X k + Xj X k Xj + A fc A| = (j^k) 

As an example of application, in analogy to Eq. ( TTTj) . we can write 

1 
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(16) 



By adding to these matrices the diagonal one 



' ei 0^ 

e 2 
1 



(17) 



we can realize the factorization 



$A 3 + B s + C 3 = AX, + B A fc + C X, 



( J VMi = l,2,3). 



(18) 



As an example, by choosing (J, k, I) = (1, 2, 3) in this equation, one obtains 
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(19) 



The matrices Ai ; 2,3 represent a generalization of Pauli matrices, but they are not the only 
possible triples. For example, we can pick the following three matrices 

( 1 \ 
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(20) 



In general, 24 triples are possible. A complete list of them can be found in Ref 



where 



an analogous problem has been treated within the context of the study of the fractional 
Pauli equation. The matrices <pj are not a generalization of Pauli matrices, as it is clear 



1 Going deeper into their algebraic structure it can be shown that they give a representation of SU(3) 
algebra generators. 
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from the computation of their commutation relations 
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We remark that the matrices in r.h.s. of Eq. ff2TT) do not provide a suitable triple for the 
cubic factorization. 

For the case of quartic root, we have found that the factorization occurs through a set of 
matrices 4x4 expressible through the Pauli matrices, and an example is given below 
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(22) 



(23) 



where 

o o y \ 10 / V 01 

We note that [xi, X2] = V%X3- We do not report all the matrices allowing the fourth order 
factorization. They can be derived by following a procedure analogous to that outlined for 
the case n — 3. 

The case of the Dirac factorization for a generic n, namely 

VA n + B n = Ar 1 + Bf 2 , (24) 

can be afforded using the same strategy as before, i. e., searching for the associated matrices 
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by means of the properties of the roots of unity. In the case of n odd we find 2 
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(25) 



where (k = 0, - ■ ■ , n — 1) are the n-th roots of unity, satisfying the conditions 



n-l 
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n-l 

J|e fc = e = 1, 

k=0 



(26) 



that determine identities of the type (fT5l) on which is based the factorization also for this 
more general case. As in the previous case the matrices Tx s 2 are not the only ones enabling 
the factorization. 

The following example, concerning the solution of a fractional evolution equation, shows 
how other matrices can naturally enter into the game. A fairly immediate consequence of 
the factorization method is that we can "linearize" the equation 



d?F(x,t)= + a n F(x, t) 



to obtain the following Dirac-like form 



dt$(x,t) = {d x Tx + aT 2 )<$>{x,t), 



(27) 



(28) 



where $ is a n-components vector. By letting $o(x) = <&(x, 0), the solution of this equation 
can be written as follows 



$(x,t)=U(t) $ (x), 



U(t) = e t{ * ld * +a * 2) . 



(29) 



Since the evolution operator U(t) is an exponential whose argument contains non-commuting 
matrices, its disentangling requires an ordered product of exponentials involving successive 
commutators of the matrices itself. At the lowest order in the commutation chain we find 3 



U(t) = e 



tfid x e ta,T2 g-(at 2 /2) [fi,f2]3aj 



(30) 



2 For n even the matrix T2 should be replaced by e % 2 n / n T2 ■ 

3 The complete disentanglement can be achieved taking into account that the matrices fj provide a repre- 
sentation of SU(N) algebra generators. 
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As discussed before (see Eq. (127|) ). the Dirac factorization method can play important 
role in the solution of functional type equations. We want to emphasize that a great deal of 
activity has been pursued in this field and that other methods can be exploited . To clarify 
this point we consider the following generalized "heat equation" 

d t F(x, t) = - Kj-dl + k F(x, t), F(x, 0) = f(x). (31) 

By taking into account the identity [3] 

POO 

exp(-cp") = de^(0exp(-c 1 />0 (p > 0, c> 0) (32) 

Jo 

where the functions <?„(£) are the so called one-sided Levy stable distribution, the solution 
of Eq. ( 131]) can be written as 

POO 

F(x,t) = / df <7i/„(0 exp [exp(t n Z%) f(x)] , (33) 

Jo 

where the part within square brackets is the kernel of the ordinary heat equation. 

In this paper we have shown how the Dirac method can be exploited to deal with frac- 
tional operators. We have also indicated the relevant importance to treat equations involving 
higher-order square root of differential operators, and provided the link with a method ex- 
ploiting the Levy distributions. We believe that the topics treated in this note are interesting 
enough to stimulate new suggestions in fractional calculus, worth to be pursued. 
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